ELSEVIER 


Available  online  at  www.sciencedirect.com 


SCIENCE 


DIRECT0 


Journal  of  Power  Sources  135  (2004)  110-121 


JOURNAL  OF 


POWER 


www.elsevier.com/locate/jpowsour 


The  effect  of  stoichiometry  on  dynamic  behavior  of  a  proton  exchange 
membrane  fuel  cell  (PEMFC)  during  load  change 

Sunhoe  Kim,  S.  Shimpalee,  J.W.  Van  Zee* 

Department  of  Chemical  Engineering,  University  of  South  Carolina,  Columbia,  SC  29208,  USA 
Received  24  February  2004;  accepted  20  March  2004 
Available  online  25  June  2004 


Abstract 

Data  are  presented  to  show  the  transient  response  of  a  proton  exchange  membrane  fuel  cell  (PEMFC)  subjected  to  change  in  the  load. 
Overshoot  and  undershoot  behaviors  of  the  steady-state  current  density  were  observed  for  various  rates  of  change  in  the  voltage  during 
constant  inlet  flow  rate  conditions.  The  results  of  these  experiments  with  a  25  cm2  triple  serpentine  flow  field  indicate  a  correlation  of 
the  overshoot/undershoot  behavior  with  initial  and  final  stoichiometry.  This  transient  analysis  is  potentially  attractive  in  the  operation  of 
vehicle  and  stationary  application  and  flow  field  designs. 
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1.  Introduction 

The  control,  design,  and  optimum  operation  of  PEMFCs 
will  require  an  understanding  of  its  transient  behavior  when 
the  current,  voltage,  or  power  changes.  These  dynamics 
would  be  important  for  residential  and  automotive  applica¬ 
tions  and  the  transient  operation  may  be  a  result  of  a  sudden 
demand  as  an  appliance  starts  or  as  a  vehicle  is  accelerated 
or  decelerated.  These  transients  may  be  of  sufficient  am¬ 
plitude  and  speed  that  fuel  flow  rates  cannot  be  adjusted 
by  feedback  control  and  thus  electrical  capacitors  are  often 
used  to  stabilize  the  system  output.  Here  we  present  exper¬ 
imental  data  to  help  understand  how  the  fuel  cell  by  default 
or  design  can  act  as  a  capacitor  during  the  power  demand 
surges. 

While  most  of  the  studies  of  PEMFCs  have  focused  on 
steady-state  behavior,  a  few  authors  have  considered  tran¬ 
sient  operation  of  stacks  [1-5],  For  example,  Hamelin  et  al. 
[1]  demonstrated  three  different  transients  in  their  work  on 
with  a  Ballard  fuel  cell  stack  model  MK5-E,  which  has  a 
total  of  35  cells,  of  10  kW  maximum  power,  connected  in 
series  with  a  cell  surface  area  of  232  cm2.  The  first  transient 
was  a  continuous  load  change  with  various  amplitudes  by 
changing  the  power  in  a  square  wave.  The  second  was  a 
higher  frequency  transient  at  very  short  times.  In  this  part 
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they  replaced  the  back-pressure  valve  with  a  metering  valve 
to  maintain  a  fixed  outlet  flow,  while  the  inlet  mass  flow 
controller  was  completely  open  to  reduce  the  transient  phe¬ 
nomena  to  the  load  alone.  Finally,  they  studied  DC  +  AC 
load  communication  and  the  noise  from  the  AC/DC  con¬ 
verter  in  the  system.  They  emphasized  the  importance  of 
studying  transient  behavior  of  fuel  cell  to  understand  how 
the  fuel  cell  stack  will  behave  under  extreme  conditions, 
for  example,  in  case  of  current  or  voltage  exceed  certain 
limits. 

Amphlett  et  al.  [2]  presented  an  analysis  of  a  hybrid 
PEMFC/battery  system  with  a  system  scaled  down  from  a 
400  kW  fuel  cell  stack  and  224  cell  lead-acid  battery  suit¬ 
able  for  submarine  systems.  They  showed  how  their  hybrid 
system  and  components  interacted  during  charging  or  dis¬ 
charging.  Those  experiments  followed  a  presentation  of  the 
dynamics  during  start-up,  current  change,  and  shut  down 
of  a  Ballard  Mark  V  35-Cell  5  kW  PEMFC  stack  where 
they  compared  the  observed  voltage  response  from  a  cur¬ 
rent  change  to  their  empirical  model  [3],  Emonts  et  al.  [4] 
studied  the  dynamics  of  a  PEMFC,  a  compact  methanol  re¬ 
former,  and  the  choice  of  a  short-term  storage  system  for  an 
automotive  application.  Kotz  et  al.  [5]  discussed  the  com¬ 
binations  of  a  PEMFC  and  supercapacitors  for  automotive 
applications.  These  experimental  studies  [1-5]  were  con¬ 
ducted  under  excess  initial  fuel  and  oxidant  conditions  and 
the  load  changes  were  not  large  enough  to  study  starved  fuel 
conditions. 
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Nomenclature 

i  current  density  (A/cm2) 

K  steady-state  gain  (A/cm2  V) 
t  time  (s) 

V  cell  voltage  (V) 

AV  cell  voltage  difference  (V) 

Greek  letters 

P  lead  time  constant  (s) 

S  time  constant  for  cell  voltage  (s) 

§  damping  coefficient 

p  lead-to-lag  ratio  (p  =  fi/x) 

x  Time  constant  (s) 


In  an  effort  to  understand  the  dynamics  of  a  PEMFC,  we 
present  data  for  a  single  cell  with  fixed  flow  rates.  The  tran¬ 
sient  operation  is  forced  by  rapid  changes  to  the  voltages 
and  this  allows  observation  of  changes  in  the  current  and 
stoichiometry  of  the  cell  as  discussed  below.  We  are  mo¬ 
tivated  by  our  recent  three-dimensional  numerical  simula¬ 
tions  for  the  transient  response  of  fuel  cell  [6,7]  in  which  we 
observed  an  overshoot  in  the  current  density  when  the  volt¬ 
age  was  changed  at  1.0  V/s.  The  dynamics  one-dimensional 
semi-empirical  model  of  Ceraolo  et  al.  [8]  also  leads  discus¬ 
sions  of  when  voltage  overshoot  observed  in  [3]  may  occur. 
Thus  we  designed  experiments  here  to  expand  an  under¬ 
standing  of  this  behavior. 


2.  Experimental 

The  objective  of  this  work  is  to  study  the  transient  cur¬ 
rent  response  when  the  cell  voltage  is  changed  at  a  rate  of 
approximately  0.2  V/s  at  fixed  feed  flow  rates.  These  fixed 
flow  rates  result  in  the  cell  being  exposed  to  different  stoi¬ 
chiometries  as  the  current  changes  as  shown  in  Table  1 .  That 
is,  here  we  present  data  for  the  response  of  the  current  when 
the  cell  voltage  was  changed  from  0.7  to  0.5  V  and  from  0.5 
to  0.7  V. 

Four  cases  can  be  described  as  changes  in  operating  con¬ 
ditions,  such  that  the  fuel  and  air  are  “normal”,  “starved”, 
or  “excess”  (see  Table  1).  For  example.  Case  la  corre¬ 
sponds  to  the  flow  rates  of  249  and  1040  cm3/min  at  standard 
conditions1  for  the  anode  and  cathode,  respectively.  These 
flow  rates  yield  a  stoichiometric  change  from  2.9/4. 8  at  r  < 
0  to  1. 2/2.0  at  t  =  oo  for  a  current  density  of  0.62  A/cm2  at 
0.7  V  and  1.48  A/cm2  at  0.5  V.  Thus  Case  la  corresponds  to 
a  change  in  operation  from  a  “excess”  to  a  “normal”  stoi¬ 
chiometry.  Case  2a  uses  the  same  flow  rates  but  we  began  the 
experiment  at  0.5  V  and  increased  the  cell  voltage  to  0.7  V 
(i.e.  0.62  A/cm2).  Case  2a  corresponds  to  a  change  in  oper¬ 


1  We  defined  the  standard  condition  as  298  K  and  101  kPa. 


ation  from  a  “normal”  to  an  “excess”  condition.  For  Case 
3a,  the  flow  rates  were  92  and  384  cm3/min,  corresponding 
to  1. 2/2.0  stoichiometry  for  the  0.7  V  initial  condition  (i.e., 
0.54  A/cm2)  and  a  stoichiometry  of  1.0/1. 7  for  t  =  oo  (i.e., 
0.5  V  and  0.63  A/cm2).  This  corresponds  to  a  change  in  the 
fuel  stoichiometry  from  a  “normal”  to  a  “starved”  condi¬ 
tion.  Note  that  for  fixed  voltage  a  single  cell  can  never  be 
completely  “starved”  but  we  use  this  as  a  description.  Prior 
to  changing  the  cell  voltage,  we  waited  30  min  to  ensure  a 
steady-state  value  of  the  current  and  to  ensure  a  well  humid¬ 
ified  membrane  at  the  voltage.  We  replicated  the  experiment 
at  least  twice  and  the  data  presented  here  demonstrates  re¬ 
producibility  of  the  transient  responses  as  discussed  below. 

All  the  data  reported  here  were  obtained  with  PRIMEA® 
Series  5510  MEA  (0.4  mg/cm2  Pt  loading,  25  p,m  nominal 
membrane  thickness,  W.L.  Gore  &  Associates,  Inc.,  Elk- 
ton,  Maryland,  USA).  The  original  active  electrode  area 
was  25  cm2,  but  sub-gaskets  for  both  anode  and  cathode 
sides  reduced  the  area  to  20  cm2.  The  gas  diffusion  me¬ 
dia  (GDM)  used  for  both  anode  and  cathode  were  16  mils 
(0.41  x  1  CP 3  m)  CARBEL™  CL.  The  cell  was  tightened 
with  8  V  at  the  torque  of  50 lb /-in. /bolt.  The  effect  of 
clamp-torque  has  been  discussed  by  Lee  et  al.  [9].  The  flow 
fields  for  both  anode  and  cathode  used  in  this  experiment 
had  triple  path  serpentine  channels  where  the  flows  were 
split  into  three  paths  at  entrance  and  converged  at  the  end. 

High  purity  hydrogen  (99.997%)  and  industrial  grade 
compressed  air  were  used.  The  fuel  cell  test  station  used 
to  control  the  electrical  load,  the  fuel  supply,  and  the  tem¬ 
perature  was  manufactured  by  Fuel  Cell  Technologies  (Los 
Alamos,  NM).  A  model  6060B  (Agilent  Technologies)  was 
used  for  the  electronic  load  bank.  Digital  mass  flow  con¬ 
trollers  (MKS  model)  were  used  to  control  the  flow  rates 
and  these  were  calibrated  with  a  bubble  flow  meter  as  dis¬ 
cussed  in  [10].  The  inlet  gases  were  bubbled  through  the 
humidity  bottles  and  from  the  temperature  of  those  bottles 
a  correlation  was  used  to  determine  the  humidity  of  inlet 
fuel  gas  [10]. 

The  operation  of  a  laboratory  scale  PEMFC  involves  inlet 
gas  and  water  vapor  diffusing  from  flow  channels  through 
the  GDM  to  electrode,  where  the  electrochemical  reaction 
takes  place.  The  non-reacted  gas  and  water  vapor  exit  the  cell 
and  pass  through  the  back-pressure  regulator  to  a  vent.  The 
back-pressure  of  both  anode  and  cathode  sides  was  101  kPa. 
The  fuel  cell  was  operated  at  the  temperature  of  70  °C  and 
the  anode  and  cathode  dew  point  temperatures  were  80  and 
70  °C,  respectively.  These  dew  point  temperatures  were  cal¬ 
culated  from  calibration  data  as  discussed  in  [10], 

The  experiment  was  conducted  in  two  parts.  First,  we 
measured  the  polarization  behavior  of  the  PEMFC  at  fixed 
stoichiometry  and  those  familiar  with  these  experiments 
will  recognize  this  indicates  that  fixed  the  flow  rates  were 
varied  according  to  the  measured  current.  These  polariza¬ 
tion  data  were  obtained  between  a  voltage  of  0.45  V  and  the 
open  circuit  voltage  in  randomized  steps  of  0.05  V.  The  flow 
rates  were  manually  changed  according  to  the  current  in  an 
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Table  1 


Flow  rates,  stoichiometries,  and  voltage  change  rates  for  experiments 


Flow 

rates,  A/C 
(cm3 /min) 

Stoichiometry 

(A/C) 

Current  density 
(A/cm2) 

Approximate 
linear  cell  voltage 
change  rate  (V/s) 

Voltage  time 
constant,  <5  (s) 

O  vershoot/undershoot(  A/cm2 ) 

t  =  0 

t  =  oo 

t  =  0 

t  =  OO 

Case  1, 
la 

0.7-0. 5  V  (excess  to  normal) 
249/1040  2.9/4.8 

1. 2/2.0 

0.62 

1.48 

0.17 

0.6 

N/A 

lb 

501/1040 

5. 8/4.8 

2. 4/2.0 

0.62 

1.48 

0.17 

0.6 

N/A 

lc 

301/1260 

3. 5/8. 7 

1. 2/3.0 

0.62 

1.80 

0.17 

0.8 

N/A 

Case  2, 
2a 

0.5-0.7V  (normal 
249/1040 

to  excess) 

1. 2/2.0 

2.9/4. 8 

1.48 

0.62 

0.17 

0.4 

N/A 

2b 

501/1040 

2. 4/2.0 

5. 8/4.8 

1.48 

0.62 

0.17 

0.4 

N/A 

2c 

301/1260 

1. 2/3.0 

3. 5/8.7 

1.80 

0.62 

0.17 

0.5 

N/A 

Case  3, 
3a 

0.7-0.5V  (normal 
92/384 

to  starved) 
1. 2/2.0 

1.0/1 .7 

0.54 

0.58 

0.22 

0.4 

0.95 

3b 

185/385 

2. 4/2.0 

1.5/1 .2 

0.57 

0.90 

0.22 

0.5 

0.95 

3c 

96/600 

1. 2/3.0 

1 .0/2.7 

0.58 

0.58 

0.22 

0.6 

1.15 

Case  4, 
4a 

0.5-0.7V  (starved 
92/384 

to  normal) 
1.0/1. 7 

1. 2/2.0 

0.58 

0.54 

0.22 

0.2 

0.39 

4b 

185/385 

1 .5/1.2 

2. 4/2.0 

0.90 

0.57 

0.22 

0.2 

N/A 

4c 

96/600 

1. 0/2.7 

1. 2/3.0 

0.58 

0.58 

0.22 

0.2 

0.35 

iterative  manner  to  maintain  the  fixed  stoichiometry.  The 
data  (Fig.  1)  were  obtained  at  three  normal  stoichiomet¬ 
ric  sets:  a  (standard),  b  (anode-rich),  and  c  (cathode-rich), 
corresponding  to  anode/cathode  stoichiometries  of  1. 2/2.0, 
2. 4/2.0,  and  1. 2/3.0,  respectively.  That  is,  for  example,  a 
stoichiometry  of  1. 2/2.0  corresponds  to  flow  rates  that  were 
1.2  times  greater  than  required  (by  the  measured  current)  for 
hydrogen  and  2.0  greater  than  that  required  on  the  cathode 
for  air. 

The  second  part  of  this  experiment  measured  the  transient 
behavior  of  the  PEM  fuel  cell.  A  two  channel  digital  oscillo¬ 
scope  (TDS  210,  Tektronix  Inc.)  was  used  to  record  simul¬ 
taneously  both  the  current  response  and  the  voltage  forcing 
function.  The  current  was  measured  using  a  “hall  effect  cur¬ 


rent  sensor”  so  that  it  would  be  converted  into  voltage  sig¬ 
nals.  Fig.  2  shows  a  schematic  diagram  of  the  experimen¬ 
tal  setup.  The  current  was  converted  into  voltage  signal  and 
then  recorded  at  the  first  channel  (CH-1)  of  the  oscilloscope. 
The  cell  potential  was  measured  directly  from  the  fuel  cell 
and  recorded  at  the  second  channel  (CH-2)  of  the  oscillo¬ 
scope.  The  cell  voltage  was  changed  after  a  steady-state  cur¬ 
rent  was  achieved  at  the  initial  voltage.  These  signals  were 
recorded  in  the  oscilloscope  from  10  s  before  the  triggering 
point,  at  t  —  0,  until  40  s  after  the  triggering  point. 

We  focused  the  transient  experiment  on  voltages  changes 
either  to  or  from  0.7  and  0.5  V.  The  test  station,  computer 
interface,  and  electronic  load  bank  limited  the  rates  of  volt¬ 
age  change  to  about  0.2  V/s.  Variation  in  this  change  was 


Fig.  1.  Steady-state  polarization  curves  for  various  stoichiometric  flows  of  the  anode/cathode:  (A)  1. 2/2.0,  (•)  2. 4/2.0,  and  (H)  1. 2/3.0. 
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Cell  Voltage 
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Fig.  2.  Schematic  diagram  of  oscilloscope  and  fuel  cell  for  transient  experiments. 


±0.05  V/s.  The  change  rates  were  not  linear  but  rather  can 
be  approximated  by  a  first-order  system.  That  is,  for  Case 
la  in  Table  1.  The  cell  voltage  change  rate  for  this  case  is 
approximately  0.17  V/s  based  on  a  linear  estimate  between 
0.7  V  at  t  —  0 s  and  0.5  V  at  t  =  1.18s.  Note,  however,  that 
the  change  is  not  linear  but  best  described  as  a  first-order 
change: 

V(f)  =  V(f  =  0)  ±  AV(1  -  e“f/5)  (1) 

where  8  is  a  time  constant  that  depends  on  the  load  char¬ 
acteristics  and  thus  on  total  current  from  the  PEMFC.  For 
Case  la,  V(t  =  0)  =  0.7  V  and  A V  =  —0.2  V  and  for  Case 
2a  V(t  —  0)  =  0.5  V  and  A  V  =  ±0.2  V.  The  time  constant, 
8,  is  sown  in  Table  1.  Also  for  Cases  3  and  4  as  shown  in 
Table  1,  the  linear  voltage  change  rate  was  approximately 
0.22  V/s  and  thus  8,  V(t)  —  0,  and  AV  follow  from  Table  1. 


3.  Results  and  discussions 

Fig.  3  shows  a  voltage-current  cycle  for  a  normal-to- 
excess-to-normal  experiment  of  Cases  la  and  2a.  Note  that 
we  performed  Case  la  followed  by  Case  2a  four  times  over 
a  period  of  approximately  4.1  h.  Fig.  4  shows  the  second 
and  third  cycles  of  Fig.  3  and  by  expanding  the  time  scale 
for  Case  la  it  allows  comparison  of  the  data  obtained  at  75 
and  135  min.  There  is  some  noise  in  the  voltage  data  corre¬ 
sponding  to  the  ±5  mV  accuracy  of  the  load  but,  in  general, 
the  experimental  data  were  reproducible.  The  response  for 
Case  la  can  be  described  as  a  first-order  (FO)  system  and 
one  does  not  observe  any  overshoot  in  the  current  density. 
That  is,  the  value  of  1.16  A/cm2  occurs  at  the  time  constant, 
r,  of  0.2  s  corresponds  to  63.2%  of  the  final  value  of  the 
current  density  difference  between  final  and  initial  values. 
That  is,  the  current  density  value  for  this  case  is 
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Fig.  3.  Overall  view  of  voltage  change  and  current  density  response  for  Cases  la  and  2a. 
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(expanded  view  with  oscilloscope),  flow  rates:  249/1040  cm3/min.  The  cell  voltage  changes  from  0.7  to  0.5  V.  Stoichiometry  change: 


i(t  =  r)  =  i(t  =  0)  +  0.632(;(f  =  oo)  -  i(t  =  0))  (2) 

where  i(t  —  0)  and  i(t  =  oo)  are  at  initial  and  final  cur¬ 
rent  density,  respectively.  The  transfer  function  of  this  FO 
response  system  to  the  voltage  change  can  be  written  as 
Laplace  domain  [11]: 

«(j)=-5tV(s)  (3) 

TS  +  1 

and  by  taking  the  inverse  Laplace  transform,  one  can  obtain 
the  time  domain  equation  for  a  step  change  in  V(s): 

i(t)  =  i(t  =  0)  +  Ki (1  -  e~,/z)(V(t)  -  V(t  =  0))  (4) 


The  gain,  K\,  is  shown  in  Table  2  to  be  negative  because 
the  forcing  function  A  V  =  V(t)  —  V(t  —  0)  is  negative  (i.e., 
0.5-0. 7  V).  Note  that  most  FO  analyses  use  a  step  function 
but  that  we  used  (Eq.  (1))  for  V(t)  from  the  experiments 
to  obtain  the  value  of  K\  with  MATLAB®  Simulink®  as 
described  in  [12].  As  discussed  below  the  response  of  current 
density  with  a  cell  voltage  change  depends  mainly  on  the 
time  constant  for  the  case  where  the  stoichiometry  is  in 
excess.  The  quality  of  the  fit  of  equation  to  the  data  is  also 
discussed  below. 

Fig.  5  shows  the  expanded  time  scale  response  for  Case 
2a,  a  change  from  0.5  to  0.7  V  at  the  same  flow  rates  as  Case 


Table  2 


Observations  and  dependent  variables  analysis  parameters 


Type  of  response 

Gain  (A/cm2V) 

Time 

constant,  z  (s) 

FO  lead  time 
constant,  f5  (s) 

Damping 
factor,  £ 

R2 

Ki  K2 

Case  1,  0.7-0. 5  V  (excess  to  normal) 

la  FO 

-4.4 

0.2 

0.9970 

lb  FO 

-4.4 

0.2 

0.9970 

lc  FO 

-5.9 

0.2 

0.9983 

Case  2,  0.5-0.7V  (normal  to  excess) 

2a  FO 

-3.6 

0.6 

0.9941 

2b  FO 

-3.6 

0.6 

0.9941 

2c  FO 

-5.1 

0.6 

0.9931 

Case  3,  0.7-0.5V  (normal  to  starved) 

3a  SO  L/L 

2.4  -2.4 

1.5 

0.4 

0.9643 

3b  FO 

-1.9 

0.02 

0.9956 

3c  SO  L/L 

2.8  -2.8 

1.1 

0.4 

0.9472 

Case  4,  0.5-0. 7  V  (starved  to  normal) 

4a  FO  L/L 

-0.1 

2.0 

19 

0.9715 

4b  FO 

-1.4 

0.4 

0.9816 

4c  FO  L/L 

-0.1 

2.0 

30 

0.9762 

The  types  of  responses  are  abbreviated:  FO,  first-order  system  response;  FO  L/L,  first-order  lead/lag  system  response;  SO  L/L,  second-order  lead/lag 
system  response  (two  gains,  K\  and  K2). 
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Fig.  5.  Case  2a  (expanded  view  with  oscilloscope),  flow  rates:  249/1040  cm3 /min.  The  cell  voltage  changes  from  0.5  to  0.7  V.  Stoichiometry  change: 
1.2/2.0-2.9/4.8. 


la.  Thus,  we  observed  a  change  from  normal  stoichiometry 
to  excess  (from  1. 2/2.0  to  2. 9/4.8  stoichiometry).  The  repro¬ 
ducibility  is  shown  to  be  good  by  rescaling  the  time  axis 
for  the  data  starting  at  45  and  105  min  and  plotting  Rep.  1 
and  Rep.  2  data  from  Fig.  3.  As  the  cell  voltage  increases, 
the  current  density  decreases  and  this  response  is  also  a  FO 
system  with  the  time  constant  corresponding  to  about  0.6  s. 
Note  that  this  FO  model  and  parameters  corresponds  to  the 
short  time  response  until  0.8  A/cm2  and  that  there  is  a  grad¬ 
ual  decrease  in  the  current  during  the  next  30  min  until  the 
value  0.62  A/cm2  is  obtained  at  0.7  V.  This  current  density 
change  for  30  min  after  a  cell  voltage  change  to  0.7  V  is  re¬ 
producible  as  shown  in  Fig.  3  and  it  is  typical  for  this  MEA 
when  there  are  transients  at  “excess”  conditions.  The  exper¬ 
iments  reported  here  do  not  yield  information  on  the  cause 


of  this  long-term  decrease.  No  undershoot  behavior  is  ob¬ 
served  for  this  case.  Case  2a  can  be  described  by  Eqs.  (3) 
and  (4)  with  the  parameters  of  Table  2.  For  Case  2a  the  cell 
voltage  change  is  a  positive  value  and  thus  the  current  den¬ 
sity  decreases  so  that  the  gain,  K\,  has  negative  value.  Again 
the  comparison  of  the  FO  model  and  the  data  is  discussed 
after  the  data  for  the  other  cases  are  presented. 

Fig.  6  shows  the  voltage  and  current  density  changes  for 
Cases  3a  and  4a.  These  results  were  obtained  over  2.0  h. 
The  current  was  unstable  (non-periodic  oscillations  were 
observed)  at  the  cell  voltage  of  0.5  V  because  the  hydro¬ 
gen  flow  approaches  starved  conditions.  This  stoichiom¬ 
etry  was  estimated  based  on  the  average  current  density 
at  0.5  V.  At  a  stoichiometry  of  1 .0  the  flow  of  hydrogen 
out  of  the  fuel  cell  will  approach  zero  because  all  of  the 
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Fig.  6.  Overall  view  of  voltage  change  and  current  density  response  for  Cases  3a  and  4a. 
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Fig.  7.  Case  3a  (expanded  view  with  oscilloscope),  flow  rates:  92/384  cm3 /min.  The  cell  voltage  changes  from  0.7  to  0.5  V.  Stoichiometry  change: 
1 .2/2.0— 1.0/1 .7. 


hydrogen  into  the  cell  will  be  reacted  before  reaching  the 
end  of  the  flow  channel.  This  is  consistence  with  our  nu¬ 
merical  predictions  [6,7]  that  indicate  for  these  conditions 
that  all  the  hydrogen  is  consumed  before  reaching  the  end 
of  the  flow  channel  and  that  not  all  the  reaction  area  is 
utilized. 

Fig.  7  expands  the  scale  of  Fig.  6  for  Case  3a.  The  sto¬ 
ichiometry  changes  from  a  normal  (1. 2/2.0)  to  a  starved 
(1.0/1. 7)  condition  (see  Table  1)  and  overshoot  behavior  is 
observed  here.  That  is,  as  the  voltage  decreases  from  0.7 
to  0.5  V,  the  current  density  increases  from  0.54  A/cm2,  at 
0.7  V,  to  about  0.95  A/cm2  after  0.22  s.  This  maximum  is 
obtained  after  the  cell  voltage  reaches  0.5  V  and  after  this 
maximum,  the  current  density  decreases  to  about  0.38  A/cm2 
for  Rep.  1  at  5.5  s  and  about  0.30  A/cm2  for  Rep.  2  at  about 
4.5  s.  Then  at  about  t  —  10  s  the  current  density  has  in¬ 
creased  from  0.38  to  0.58  A/cm2  for  Rep.  1.  Rep.  2  shows 
a  slightly  different  behavior  where  the  current  density  in¬ 
creases  from  0.30  A/cm2  to  an  apparent  second  overshoot 
of  0.68  A/cm2  at  about  t  =  10  s.  Typically,  we  could  not 
confirm  that  this  apparent  second  overshoot  was  different 
from  noise  but  this  second  overshoot,  shown  in  Fig.  7,  cor¬ 
responds  to  the  maximum  value  we  observed  over  many 
replicates.  Fig.  7  also  shows  second  undershoot  behavior  at 
about  t  =  17  s  which  again  was  the  extreme  that  we  ob¬ 
served  over  many  replicates.  Note  that  there  is  some  un¬ 
dershoot  with  the  cell  voltage  in  the  first  3  s  but  that  this 
undershoot  was  only  slightly  responsible  for  the  width  of 
the  overshoot  current  peak  as  determined  by  other  experi¬ 
ments. 

We  now  explain  the  overshoot  behavior:  Since  the  elec¬ 
trochemical  reaction  is  greater  at  0.5  than  at  0.7  V  and  be¬ 
cause  the  potential  of  the  electrodes  follows  the  cell  voltage 
without  any  appreciable  time  lag,  the  local  reaction  (i.e., 
current)  uses  the  excess  hydrogen  and  oxygen  without  limi¬ 


tation.  (The  time  constants  for  limitations  are  discussed  be¬ 
low.)  This  unlimited  use  continues  until  0.95  A/cm2  (based 
on  20  cm2  of  MEA  area).  The  current  density  starts  to  de¬ 
crease  at  0.22  s  because  all  of  the  “excess”  hydrogen  in  the 
GDM  and  flow  channels  is  consumed.  Note  that  hydrogen  is 
still  flowing  into  the  cell  but  that  the  distribution  of  current 
is  highly  non-uniform  at  0.22  s  and  becomes  more  uniform 
as  discussed  in  [7].  Note  also  that  the  overshoot  is  greater 
than  the  20%  that  would  be  expected  based  on  the  average 
stoichiometry.  One  might  be  concerned  that  the  MEA  hy¬ 
dration  may  also  limit  reaction  but  our  experience  and  mea¬ 
surements  indicate  that  dehydration  occurs  at  a  slower  time 
scale  than  the  response  shown  in  Fig.  7  for  this  MEA. 

The  undershoot  behavior  for  Case  3a  can  be  explained  by 
considering  the  current  distribution.  That  is,  we  hypothesize 
that  the  overshoot  behavior  results  in  a  highly  non-uniform 
current  distribution  so  that  the  consumption  of  hydrogen  al¬ 
lows  for  a  significant  portion  of  the  PEMFC  on  the  hydro¬ 
gen  side  to  contain  ambient  air.  Ambient  air  flows  into  the 
hydrogen  side  to  equilibrate  the  atmospheric  pressure  much 
like  a  spring  in  a  spring-dashpot.  We  label  this  a  “vacuum” 
effect.  This  description  requires  the  highly  non-uniform  cur¬ 
rent  distribution  at  the  peak  current  density  to  become  more 
uniform  after  10  s  as  the  current  density  decreases  from  the 
peak  and  it  requires  time  to  expel  the  ambient  air  and  use 
all  of  the  electrode  area.  The  time  between  the  crossing  un¬ 
der  the  final  current  density  values  and  the  return  to  the  final 
current  density  is  approximately  9  s  for  Rep.  1  and  7  s  for 
Rep.  2.  The  Case  3a  can  be  described  with  a  second-order 
lead/lag  (SO  L/L)  current  response  model  as  discussed  be¬ 
low. 

Fig.  8  shows  the  expanded  time  scale  behavior  for  Case  4a, 
a  change  from  0.5  to  0.7  V  at  the  same  flow  rates  as  Case  3a. 
Thus  we  observe  the  behavior  as  the  stoichiometry  changes 
from  a  starved  to  normal  condition  (i.e.,  stoichiometry: 
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Fig.  8.  Case  4a  (expanded  view  with  oscilloscope),  flow  rates:  92/384  cm3 /min.  The  cell  voltage  changes  from  0.5  to  0.7  V.  Stoichiometry  change: 
1.0/1. 7-1. 2/2.0. 


1.0/1 .7  to  1. 2/2.0).  Here,  we  observe  undershoot  behavior 
in  the  current  density.  At  the  cell  voltage  of  0.5  V,  which  is 
a  staved  condition,  the  reaction  is  probably  not  uniform  and 
the  active  area  utilization  is  not  complete.  However,  when 
the  cell  voltage  is  changed  to  0.7  V,  the  electrochemical 
driving  force  is  less  and  the  current  density  decreases.  The 
lower  current  density  yields  unused  fuel.  Then,  the  unused 
fuel  expels  the  ambient  air  and  the  non-utilized  reaction 
area  becomes  exposed  to  hydrogen.  This  re -exposure  leads 
to  an  increase  of  current  density  after  about  1  s.  The  time 
for  the  recovery  of  the  undershoot  behavior  is  about  8  s. 
The  responses  of  Case  4a  can  be  described  with  a  first-order 
lead/lag  (FO  L/L)  current  response  model  as  describe  below. 

Fig.  9  shows  a  comparison  of  the  Cases  of  la-lc  and  it 
allows  analysis  of  the  stoichiometry  effect.  The  three  cases, 


a,  b,  and  c,  show  similar  first-order  trends.  The  response 
for  Cases  la-lc  can  be  classified  as  a  first-order  system 
because  one  does  not  observe  an  overshoot  in  the  current 
density  when  the  fuel  is  excess.  The  final  performance  of 
Case  lc  (cathode-rich  stoichiometry)  shows  the  highest  cur¬ 
rent  density  (1.80  A/cm2)  and  the  initial  current  densities 
are  equal  indicating  ohmic  limitations  at  the  high  initial  sto¬ 
ichiometries  consistent  with  Fig.  1  shown  in  Table  1.  Cases 
la  and  lb  show  identical  current  density  responses  (0.65  and 
1.48  A/cm2  for  initial  and  final  values,  respectively).  The 
linear-approximated  cell  voltage  change  rates  for  these  three 
cases  are  0.17  V/s  and  the  time  constant,  8 ,  for  Eq.  (1)  are 
close  for  all  three  cases.  The  FO  model  for  the  response  of 
the  current  and  a  FO  response  occur  for  all  of  these  cases 
because  they  are  not  hydrogen-staved  condition.  Table  2 


Time  (s) 


> 

(D 

Ul 

(0 

+* 

o 

> 

a> 

O 


Fig.  9.  Comparison  of  the  stoichiometry  effect  for  Cases  la-lc  (excess  to  normal).  See  Table  1  for  flow  rates. 
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Fig.  10.  Comparison  of  the  stoichiometry  effect  for  Cases  2a-2c  (normal  to  excess).  See  Table  1  for  flow  rates. 


shows  that  the  time  constant  for  each  case  is  the  same  and 
that  the  gain  is  controlled  by  the  air  stoichiometry.  This  gain 
could  be  obtained  from  steady-state  polarization  curves  such 
as  those  shown  in  Fig.  1  and  thus  for  this  configuration  of 
MEA,  GDL,  gasket,  flow  field,  and  membrane  water  content 
the  FO  response  can  be  calculated  for  design  purposes.  Fu¬ 
ture  work  will  quantify  the  effect  of  A  V  on  r  and  K\  for  the 
excess  stoichiometry  and  other  configurations.  The  lack  of 
a  starved  condition  conditions  yields  no  overshoot  behavior. 
The  R 2  shown  in  Table  2  is  one  measure  of  good  agreement 
between  the  FO  model  and  the  data  and  Ref.  [12]  contains 
graphs  of  these  comparisons. 

Fig.  10  shows  a  comparison  of  Cases  2a-2c.  These  graphs 
show  similar  FO  current  density  responses  where  the  current 
density  decreases  when  the  cell  voltage  increases.  The  initial 
performance  of  the  Case  2c  is  higher  than  that  of  other  two 
cases  because  the  higher  cathode  stoichiometry  air  reduces 
mass  transfer  limitations.  This  response  system  can  be  char¬ 
acterized  with  FO  response  of  exponential  decay  because 
there  is  excess  fuel.  The  FO  model  with  the  parameters  in 
Table  2  agrees  with  the  experimental  data  and  no  undershot 
behavior  is  observed.  Note  that  although  the  time  constants 
are  equal,  the  gain  reflects  the  initial  stoichiometry.  Note 
also  that  these  parameters  correspond  to  short  time  since 
the  final  approach  to  steady-state  requires  about  30  min  (see 
Fig.  3  and  discussion  above).  This  slow  long-term  decay¬ 
ing  is  responsible  for  the  difference  in  K\  for  the  respective 
conditions  of  Cases  1  and  2.  Thus  K i  =  —5.9  A/cm2  V  for 
Case  lc  but  K\  —  —5.1  A/cm2  V  for  Case  2c. 

Fig.  11a  shows  a  comparison  of  Cases  3a-3c.  These  data 
have  more  noise  due  to  the  “starved”  condition  but  the  data 
are  sufficiently  reproducible  to  compare.  Case  3b  does  not 
show  the  overshoot/undershoot  behavior  because  there  is  ex¬ 
cess  hydrogen.  Cases  3a  and  3b  start  at  the  same  current 
density  at  0.7  V  (consistence  with  Fig.  3)  but  only  Case  3b 
shows  a  FO  responses  as  indicated  in  Table  2.  Note  that  the 


small  time  constant  of  0.02  s  indicating  the  lack  of  limita¬ 
tions  consistent  with  the  lack  of  mass  transfer  or  membrane 
hydration  issues.  Case  3c  has  the  higher  air  stoichiometry 
than  Case  3a  and  thus  the  initial  current  density  for  Case 
3c  is  slightly  the  larger  (0.58  versus  0.54  A/cm2).  The  re¬ 
sponses  of  Cases  3a  and  3c  can  be  classified  as  a  modified 
second-order  response  called  as  second-order  lead/lag  sys¬ 
tem  (SO  L/L).  The  details  of  this  modification  can  be  found 
in  [12].  The  model  can  be  written  in  the  Laplace  domain: 

'M  =  (^+lV+K2U+i)VW  <5) 

The  time  domain  function  for  this  system  can  be  described 
by  the  sum  of  two  current  densities,  i\  and  A,  where  4  cor¬ 
responds  to  the  current  density  for  all  times  beyond  the  time 
where  the  current  density  is  a  maximum  and  is  described: 

r2^  +  2^  +  i,  =  K(V(t)  -  V(t  =  0))  (6a) 

and  where  i2  corresponds  to  the  current  density  at  times 
before  and  at  the  maximum  and  is  described: 

i2  =  K2  [V(t)  -  V(t  =  0)]  (6b) 

MATLAB®  [13]  was  used  to  fit  the  Eq.  (5)  to  the  data  and 
obtain  the  parameters.  The  parameters,  steady-state  gains, 
K\  and  K2,  the  time  constant,  r,  and  the  damping  factor, 
§,  listed  in  Table  2.  When  the  cell  voltage  is  decreased  to 
0.5  V,  AV  =  —0.2  V,  and  the  current  density  increases  in¬ 
stantaneously  to  0.95  and  1.15  A/cm2  for  Cases  3a  and  3c, 
respectively.  We  model  this  instantaneous  increase  with  the 
gain  K2.  The  condition  of  K2  with  AV  is  explored  in  our 
next  paper  [12].  The  cell  voltage  difference,  AV  has  nega¬ 
tive  value  in  this  case,  —0.2  V,  and  K2  has  negative  value  so 
that  the  current  density  response  has  positive  value.  After 
the  maximum  current  density,  all  the  hydrogen  stored  in  the 
cell  is  used  and  the  current  density  begins  to  decrease.  This 
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Fig.  11.  (a)  Comparison  of  the  stoichiometry  effect  Cases  3a-3c  (normal  to  starved).  See  Table  1  for  flow  rates,  (b)  Fitting  of  second-order  lead/lag 
response  for  Case  3a.  (c)  Fitting  of  second-order  lead/lag  response  for  Case  3c. 


decrease  is  modeled  with  the  second-order  part  of  Eq.  (5) 
and  (6).  The  gain,  K\ ,  has  positive  value  and  the  AV  has 
negative  value  so  that  the  current  begin  to  decay  from  the 
overshoot  peak.  The  time  constant  can  be  obtained  at  63.2% 
of  this  decrease.  The  current  density  value  at  t  =  0  in  that 
equation  in  this  case  is  the  overshoot  peak  and  that  of  t  = 
oo.  The  undershoot  behavior  is  affected  by  the  damping  co¬ 
efficient  in  the  second-order  term.  The  time  constants,  r,  for 
Cases  3a  and  3c  are  1.5  and  1.1s,  respectively.  This  indicates 
that  the  hydrogen  consumption  is  about  1.4  times  faster  for 
Case  3c.  These  parameters  could  be  used  for  interpolation 


for  any  cathode  stoichiometry  ratio  between  these  two  cases. 
The  instantaneous  gains,  AC,  for  Cases  3a  and  3c  are  dif¬ 
ferent  based  on  the  difference  in  overshoot  at  t  — »■  0  listed 
in  Table  2.  Note  that  the  R2  values  are  somewhat  smaller 
for  Cases  3a  and  3c  indicating  less  agreement  between  the 
equation  and  the  data.  Fig.  lib  and  c  shows  that  these  com¬ 
parisons  are  adequate  for  a  description  of  the  dynamics. 

It  is  interesting  to  note  that  the  overshoot  peak  at  about 
0.95  A/cm2  for  Case  3a  is  only  slightly  larger  than  the  value 
of  about  0.90  A/cm2  observed  for  Case  3b.  The  hydrogen  is 
not  starved  for  Case  3b  (stoichiometry  change  from  2.4  to 
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1.5  for  the  anode)  and  thus  one  might  believe  that  the  max¬ 
imum  peak  is  related  to  the  ohmic  drop  across  the  MEA. 
Thus  Case  3b  is  controlled  by  oxygen  electrode  and  the  IR 
since  there  is  sufficient  hydrogen.  Thus  Case  3b  gives  the 
stable  cell  performance  at  0.5  V,  while  other  two  cases  show 
“vacuum”  effects  and  unstable  oscillations  due  to  insuffi¬ 
cient  hydrogen.  Also,  the  undershoot  behavior  due  to  the 
“vacuum”  effect  after  overshoot  peak  for  Case  3c  is  more 
significant  than  that  of  Case  3a,  due  to  the  cathode-rich  sto¬ 
ichiometry.  The  overshoot  peak  values  for  Cases  3a  and  3c 
are  about  0.95  and  1.15  A/cm2,  respectively.  The  undershoot 
behaviors  were  shown  at  about  t  =  5  s  with  a  depth  of  about 
0.37  A/cm2  for  Case  3a  and  at  about  t  =  2.5  s  and  0.2  A/cm2 
of  undershoot  value  for  Case  3c.  With  similar  hydrogen  flow 
rates  for  Cases  3a  and  3c  are  similar,  92  and  96  cm3/min,  re¬ 
spectively,  the  hydrogen  is  consumed  faster  for  Case  3c  than 
for  Case  3a  because  the  oxygen  is  not  limited.  The  electro¬ 
chemical  reaction  is  more  limited  by  oxygen  mass  transfer  at 
cell  voltage  of  0.5  V  than  0.7  V.  Abundant  oxygen  for  Case 
3c  results  in  faster  consumption  of  hydrogen  than  Case  3a 
and  a  larger  “vacuum”  effect  with  more  ambient  air  flow¬ 
ing  in.  Also  the  second  overshoot  is  observed  in  Case  3c 
because  we  used  Rep.  2  of  Fig.  7  in  Fig.  11a.  Again  while 
we  cannot  completely  distinguish,  the  second  peak  observed 
at  about  t  —  10  s  with  the  value  of  about  0.62  A/cm2  from 
the  noise,  one  would  expect  that  a  second  oscillation  cor¬ 
responding  to  the  dashpot  system  mentioned  above  is  more 
significant  with  excess  air. 

Fig.  12  shows  a  comparison  of  the  Cases  4a-4c.  Case 
4b  does  not  show  undershoot  behavior  because  the  hydro¬ 
gen  flow  is  not  starved.  Also,  the  initial  current  density, 
0.90  A/cm2,  is  the  same  as  the  final  of  Case  3b  from  Fig.  11a. 
The  transient  behavior  of  Case  4b  can  be  classified  as  a 
FO  response  system  as  explained  above.  The  final  current 
density  value  in  Fig.  12  is  about  0.6  A/cm2  for  all  cases. 


However,  the  undershoot  behavior  of  cathode-rich  condition, 
Case  4c,  shows  similar  behavior  as  that  shown  in  Case  4a. 
The  undershoot  depth  for  Cases  4a  and  4c  are  about  0.40 
and  0.38  A/cm2  and  the  recovery  times  for  both  cases  are 
also  very  similar  with  similar  hydrogen  flow  rates.  We  ex¬ 
plain  this  similarity  as,  the  cell  performance  at  0.7  V  is  more 
affected  by  ohmic  limitations  rather  than  oxygen  mass  trans¬ 
fer  limitation.  Thus  the  extent  of  recovering  from  hydrogen 
starved  conditions  is  similar  for  both  cases.  This  similarity 
strongly  supports  the  hypothesis  that  the  undershoot  behav¬ 
ior  and  recovery  rate  depends  on  the  hydrogen  utilization. 

Cases  4a  and  4c  can  be  described  by  a  first-order  lead/lag 
(FO  L/L)  response  of  current  density.  The  FO  L/L  system 
can  be  written  in  the  Laplace  domain  as: 

i(s)  =  K\  ^2  Vis)  (7) 

rs  +  1 

The  time  domain  of  the  transfer  function  above  is 

i(t)  =  i{t  =  0)  +  KX{  1  -  (1  -  p)  e-?/T)(R(0  -  Vit  =  0)) 

(8) 

where  p  is  the  lead-to-lag  ratio  (—P/r)  [11,13].  The  pa¬ 
rameters,  lead  time  constant,  ft,  lag  time  constant,  r,  and 
steady-state  gain,  K\ ,  are  listed  in  Table  2.  The  gain,  K \ ,  and 
lead  time  constant,  ft,  affect  the  instantaneous  undershoot 
response.  The  gain  represents  the  difference  between  initial 
and  final  values.  Although  the  initial  current  density  values 
for  Cases  4a  and  4c  shown  in  Table  1  are  0.58,  the  aver¬ 
age  current  density  at  0.5  V  as  mentioned  in  Fig.  8,  for  both 
cases,  the  model  fittings  were  calculated  based  on  Fig.  12. 
The  undershoot  peak  is  determined  by  lead-to-lag  ratio,  p 
and  K\ .  The  A  V  has  positive  value,  0.2  V,  and  the  gain  has 
negative  value  so  that  the  current  response  is  negative  with 
the  cell  voltage  change.  Then  the  lag  time  constant  in  this 
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Fig.  12.  Comparison  of  the  stoichiometry  effect  in  Cases  4a-4c  (starved  to  normal).  See  Table  1  for  flow  rates. 
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case,  r,  affects  the  current  density  makeup.  The  lag  time  con¬ 
stants,  r,  were  obtained  from  experimental  data  by  applying 
Eq.  (2).  The  initial  current  density  in  this  case  is  the  under¬ 
shoot  peak  value.  While  the  other  two  cases  show  first-order 
lead/lag  system.  Case  3b  shows  a  modified  first-order  expo¬ 
nential  decay.  The  hydrogen  stoichiometry  is  changed  from 
1.5  to  2.4  for  Case  3b. 

4.  Conclusions 

Data  were  presented  to  show  the  current  density  response 
when  the  cell  voltage  was  changed  rapidly  for  fixed  flow 
rates.  The  responses  depend  on  the  relative  stoichiometry 
at  t  =  0  and  oo.  Overshoot,  undershoot,  and  overshoot/ 
undershoot  behavior  was  observed.  For  excess  stoichiom¬ 
etry,  the  response  was  modeled  with  a  FO  equation.  Thus 
no  overshoot  or  undershoot  behaviors  were  observed 
for  the  excess  stoichiometries.  For  “starved”  conditions, 
pseudo-second-order  behavior  was  observed  and  classified 
with  lead/lag  models.  Parameters  for  these  models  were 
presented  and  the  models  agree  well  with  the  experimental 
data.  Undershoot  behaviors  were  observed  when  hydrogen 
stoichiometry  is  changed  from  starved  to  normal  condi¬ 
tion.  The  undershoot  behavior  was  explained  in  terms  of  a 
non-uniform  hydrogen  controlled  current  distribution.  The 
undershoot  peak  is  affected  by  the  gain,  K,  and  lead-time 
constant,  /3,  with  the  cathode  stoichiometry.  Overshoot 
behavior  was  explained  in  terms  of  excess  hydrogen  con¬ 
sumption,  non-uniform  hydrogen  controlled  current  density, 
and  a  sufficient  consumption  of  hydrogen  to  draw  ambient 
air  into  the  flow  field  at  the  exit  of  the  cell.  These  over¬ 
shoot  peaks  are  affected  by  the  gain,  Ki,  for  instant,  current 
density  increase  with  cell  voltage  decrease. 

Finally,  overshoot  followed  by  undershoot  behavior  was 
observed  and  explained  in  terms  of  hydrogen  replacing  the 
ambient  air  as  the  current  density  distribution  became  more 
uniform.  The  undershoot  following  the  overshoot  peak  is  af¬ 
fected  by  second-order  gain,  K\,  for  the  steady-state  value 
at  t  =  oo  and  damping  coefficient,  if  for  its  magnitude. 
The  cathode  stoichiometry  affected  the  time  constant,  r,  and 


magnitude  of  overshoot  peak.  It  supports  that  the  overshoot 
or  undershoot  behaviors  were  mainly  dependent  on  hydro¬ 
gen  utilization. 
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